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Abstract:  By  expressing  the  transient  electromagnetic  behaviors  in  terms  of  continuous  orthonormal  Laguerre 
polynomials,  we  can  obtain  an  unconditionally  stable  solution  of  the  time  domain  electric  field  integral  equation 
(TD-EFIE)  for  three-dimensional  (3-D)  arbitrary  shaped  conducting  bodies.  Besides  using  Gaussian  type  pulses, 
rectangular  and  triangular  pulses  are  also  used  as  incident  waves  in  this  method.  However,  because  of  the 
discontinuity  in  a  rectangular  pulse,  Gibbs  phenomenon  will  occur  around  the  point  of  discontinuity  when  a 
continuous  basis  functions  are  used  to  approximate  the  incident  wave  in  a  least  square  sense.  Noting  that  we 
deal  with  discrete  data  during  our  computation  in  a  computer,  we  introduce  the  discrete  Laugerre  functions  to 
solve  TD-EFIE.  They  are  exactly  othonormal  in  a  discrete  sense.  In  this  paper,  we  use  the  discrete  Laguerre 
basis  functions  to  approximate  its  continuous  counterparts  and  then  use  them  to  express  the  rectangular  incident 
wave  and  the  response.  Simulation  results  show  that  there’s  no  Gibbs  phenomenon.  Furthermore,  the 
computation  of  the  Laguerre  transform  of  the  incident  wave  is  more  efficient. 
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1.  TD-EFIE 


Let  S  denote  the  surface  of  an  electric  conducting  body  illuminated  by  a  transient  electromagnetic  wave. 
Since  the  total  tangential  electric  field  is  zero  on  the  conducting  surface  for  all  times,  we  have 


— A(r,  t)  +  V0( r,  t) 
o  t 


=  [Ei(r,0L.reS 


(1) 


where  E'  is  the  incident  field  and  A  and  0  are  the  magnetic  vector  potential  and  the  electric  scalar  potential 
given  by 


J(r',r)rfy 

(2) 

R 

9(r',r)rfs, 

(3) 

R 

In  (2)  and  (3),  R  —  |  r  —  r'  |  represents  the  distance  between  the  arbitrarily  located  observation  point  r  and  the 
source  point  r'.  r  =  t  —  R  /  c  is  the  retarded  time,  q  is  the  surface  charge  density.  //  and  £  are  the  permeability 
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and  the  permittivity  of  the  medium,  and  c  is  the  velocity  of  the  propagation  of  the  electromagnetic  wave  in  that 
space.  Equation  (1)  with  (2)  and  (3)  constitutes  a  TD-EFIE  from  which  the  unknown  current  J  may  be 
determined. 

The  surface  of  the  structure  to  be  analyzed  is  approximated  by  planar  triangular  patches.  We  define  the 
spatial  basis  function  associated  with  the  n -th  common  edge  as 

+  [-'■ 

f«  (r)  =  C  (r)  +  f„  (r)  ,  f*(r)  = 

[  0,  r  eT* 

where  ln  and  A*  are  the  length  of  the  edge  and  the  area  of  triangle  T* .  p  *  is  the  position  vector  defined  with 
respect  to  the  free  vertex  of  T* .  The  electric  current  J  on  the  scattering  structure  may  be  approximated  in 
terms  of  the  vector  basis  function  as 


— — n  *  t  gT± 

2  A~"’ 


(4) 


J(r,0  =  £j„(0f„(r) 


(5) 


where  N  represents  the  number  of  common  edges,  discounting  the  boundary  edges  in  the  triangulated  model  of 
the  conducting  object.  We  introduce  a  new  source  vector  e(r,  t)  defined  by 


By  using  (5)  and  (6),  we  may  express 


J(r,  0  =  T-e(r,  t) . 

6 1 


e(r,  t)  =  Yjen(t)  f„(r). 


(6) 


(7) 


We  now  solve  (1)  by  applying  Galerkin’s  method  in  the  MoM  context  and  hence  the  testing  functions  are  same 
as  the  expansion  functions.  By  choosing  the  spatial  expansion  function  fm  (r)  also  as  the  spatial  testing 
functions,  we  have  from  (1) 


<  f,„  (r),  —  A(r,  0>  +  <  f„,  (r),  V<P(r,  t)  >=<  f,„  (r),  E'  (r,  t)  > 
o  t 


(8) 


where  m  =  1,  N  .  In  computing  (8),  we  assume  that  the  unknown  transient  quantity  does  not 
change  appreciably  within  the  triangle  so  that 

R  Rpq  i  i 

T  =  t~—  ->  Tpq=t--^~  Rpq  =  \rcp  -rcq\  19) 

L  1  7  L  mn  1  ’  lvmn  1  m  ln  \y ) 

C  C  1  1 

where  p  and  q  are  +  or  -.  rjf  is  the  position  vector  of  the  center  in  triangle  T* .  With  the  assumption  (9),  (8) 
can  be  written  as 


II 

«= i  ps  L 


d 2  bpq 

uapq  —^e  (rpq)  +  ^-e  (rpq) 
at  s 


=  vEm(t) 


where 


apq  = 

mn 


An 


W) 


dS'dS . 


bpq  =  — 
An 


dS'dS . 


VE(t)=[fm(r)-E\r,t)dS. 

•o 


(10) 

(11) 

(12) 

(13) 
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2.  Continuous  Laguerre  Polynomials  and  the  formulation 


An  orthonormal  basis  function  set  can  be  derived  from  the  Laguerre  functions  through  the  representation  ([2], 

[3]) 

</>j(st)  =  ^e~s,,2L/.(st)  (14) 

where  Lj{t)  is  the  Laguerre  polynomial  of  degree  j  and  s  is  the  scaling  factor.  The  continuous  Laugerre 
transform  for  e„(t)  can  be  given  as 


oo  oo 

en,j  =  jen(t)</>j(st)dt ,  e„(t)  =  (15) 

o  j=o 

The  expressions  of  expanding  the  first  and  the  second  derivative  of  the  transient  coefficient  are  given  as, 
respectively, 


d2 

dt2 


d  ^ 

oo 


2  "J 

1 

4  tn’J 


4M) 


i  j 

V  +  IX* 
A=0 

7-1 

j= 0  L  "t  i=0 


<pXst). 


(16) 

(17) 


Substituting  (17)  and  (15)  into  (10)  and  taking  a  temporal  testing  with  (j),{st )  ,  we  have 


III 

«=1  PA  7=0 


C s1  bpq A 

—  uapq+  — 

4r“~mn 

s  j 


en,j+S 


lVaP:Jt(j~k)en,k 

k= 0 


I; 


(  npq  \ 

^ mn  I  _  y-E 

m,i 


where 


Rpq 

,  mn 


J*00 

0 


rm 

3t-s-^-\d(st) 


(18) 


(19) 


VE  =  \jXst)V°(t)d(st).  (20) 

We  note  that  I  j  —  0  when  j  >i .  Therefore  we  can  write  the  upper  limit  of  the  third  summation  symbol  as  i 
instead  of  oo  in  (18).  In  this  result,  moving  the  terms  including  en  . ,  which  is  known  for  j  <i ,  to  the  right- 
hand  side,  we  obtain 


zz 

n= 1  p,q 


(  1  UPq^ 

—RaZ  +  — 

V4  £  j 


e„  J, 


Ri 


V  c  J 


N  (  s 2 


=  Ki-  ZZZ 

«=1  p,q  y=0 


bpq 

—paZ+  — 

V4  *  y 

7-1 


^  —El 

v  c  y 


■ZZZ's'2/<«  Z  (7  - 

n= 1  p,q  7=0  A:=0 


pq  A 

mn 


Rewriting  (21)  in  a  simple  form,  we  have 


YaEe  -  VE  +PE 

r  m,z  m,  z 


(21) 


(22) 


where 
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a„ 


=  Z 


f  s2  bpq 

—  Liam+  — 

4r*  mn 
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\ 

exp 

_ ^  vmn 

/ 

V  2c  7 

pe. =-y  y 
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1 

z 
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7? 


(23) 


(24) 


It  is  important  to  note  that  [«fn  ]  is  not  a  function  of  the  degree  of  the  temporal  testing  function.  Therefore,  we 


can  obtain  the  unknown  coefficients  by  solving  (22)  through  an  increase  in  the  degree  of  the  temporal  testing 
functions. 


3.  The  discrete  Laugerre  basis  functions  and  the  transforms  for  rectangular  pulse 


The  discrete  Laguerre  functions,  defined  in  the  Z  domain,  can  be  written  as  ([1]) 

'¥J{z,a)  =  J^a2  (Q~Z_  ;  |a|<l,  j  =  0,1,2,.. 


(25) 


(1  -az~l)J+l 

The  constraint  on  the  pole  |a|<l  is  set  to  make  the  functions  causal  and  stable.  It  has  the  following  recursive 
form  in  time  domain 

¥j+ 1  (k)  =  a  y/j  ( k )  -  y/ y  (k  - 1)  +  a  y/ j+x  ( k  - 1)  (26) 

The  discrete  Laguerre  functions  are  orthonormal  in  a  discrete  sense,  which  means 


oo  i  J 

Tj  -  X  V,  ( k ,  d)y/j  (k,  a)  =  —  ^  'P,  (z,  a)^  (1  /  z,  a )  —  =  Sj 
k= o  2^7  z 


(27) 


The  integration  is  along  the  unit  circle  in  Z  plane.  In  the  limit  At  approaches  zero,  The  sample  value 
<f>j(sk  ■  At)  and  y/ j(k,a)  become  equal  as  they  are  related  by 


a  =  e  2  ,  or  5  = - lna 


1_ 

At 


(28) 


where  At  is  the  sampling  interval.  The  first  5  orders  of  continuous  and  discrete  Laguerre  functions  for  a  =  0.9 
and  At  =  1  are  shown  in  Fig.  1,  where  the  lines  are  the  continuous  functions  and  dots  are  the  discrete  ones. 

The  discrete  Laguerre  transform  can  be  defined  as 


dnj  =  z  e»  (*■  Vy  (*» «) .  (k)  =  Z  d»jV/j  (k> a ) 

k=Q  7=0 

where  e„(k)  =  e„(k  At)  is  the  Mi  sample  of  e„(t). 

Using  the  discrete  Laguerre  basis  functions  to  approximate  the  continuous  ones,  (20)  becomes 

OO 

C=£  ^(kWmE(k) 


(29) 


(30) 


The  integration  is  replaced  by  the  summation.  This  approach  can  reduce  the  computational  load  and  can  avoid 
errors  introduced  in  numerical  integration. 

The  continuous  Laugerre  transform  of  a  rectangular  pulse  can  be  analytically  computed  as 

2  i 

rJ{st)dt  = - <pj{sm  -2 y  q. 


C7  = 


(31) 


k= 0 


where  t\  and  h  are  the  start  and  stop  time  for  the  pulse  with  an  amplitude  of  IV.  The  discrete  Laguerre 
transform  of  a  rectangular  pulse  can  be  computed  as 
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d  _  ^  (/:)  _ aWki) - rtfi)]- k-i(^i) ~ Vj-  1^2)]  d 


1  -  a 


j- 1 


(32) 


where  k\  and  k2  are  the  sampled  values  at  the  start  and  stop  time.  As  a  reference,  we  compare  the  coefficients  of 
the  Laguerre  transforms  with  that  from  a  Fourier  series.  It’s  known  the  Fourier  series  for  a  rectangular  pulse  can 
be  written  as 


_2(t2-tl)  _ -sin(2k7rtl/T)  +  sin(2k7rt2/T)  _cos(2kKtlIT)-cos(2kKt2IT) 

0  T  ’  J  kn  ’  J  kn 

The  transform  data  for  the  three  kinds  for  a  rectangular  pulse  with  t\  =  1  lm  (light  meter),  t2  =  2  lm,  time  span  T 
=  10 lm  are  shown  in  Fig.  2.  We  can  see  that  for  large  orders,  the  values  of  the  discrete  Laugerre  transform 
approaches  to  zero  exponentially  while  the  other  2  are  oscillating. 

Reconstruction  of  the  original  signal  by  a  discrete  Laguerre  transforms  is  more  convergent  than  the  continuous 
transforms,  which  is  also  shown  in  Fig.  3.  Unlike  the  other  2  transforms,  the  reconstructed  signal  by  the 
discrete  transform  doesn’t  display  the  well-known  Gibbs  phenomenon. 


4.  Numerical  Simulation 

In  this  paper,  an  example  of  scattering  analysis  of  a  dipole  is  shown  as  in  Fig.  4.  The  dipole  is  a  conducting 
strip  with  one  meter  long.  A  rectangular  incident  plane  wave  in  the  time  domain  as  shown  in  Fig.  3  is 
propagating  along  the  X  axis.  We  analyze  the  induced  current  at  the  middle  edge  of  the  dipole.  Figure  5  shows 
the  comparison  of  the  currents  obtained  by  continuous  and  discrete  Laugerre  basis  functions.  When  we  enlarge 
the  results  in  early  time  in  Fig.  6,  we  can  see  the  response  obtained  by  discrete  Laugerre  functions  is  zero  when 
there’s  no  input  incident  wave.  It’s  more  reasonable  than  that  obtained  by  the  continuous  basis  functions.  The 
oscillation  shown  in  the  early  time  for  the  continuous  basis  functions  is  caused  by  the  Gibbs  phenomenon. 

5.  Conclusions 

In  this  paper  we  present  the  discrete  Laugerre  basis  functions  in  the  solution  of  TD-EFIE.  The  discontinuous 
rectangular  incident  wave  is  used  in  the  scattering  analysis.  The  discrete  basis  functions  are  superior  to  its 
continuous  counterpart  in  the  fact  of  representation  the  incident  wave.  It’s  more  efficient  in  computing  the 
transform  data  and  the  Gibbs  phenomenon  doesn’t  exists  if  we  use  the  discrete  basis  functions.  Since  the 
discrete  Laguerre  functions  cannot  approximate  its  continuous  counterpart  in  higher  orders,  more  work  is  to  be 
done  for  the  improvement  of  the  algorithm. 
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Discrete  and  continuous  Laguerre  functions  for  a=0.9 


coeficients  for  a=0.59652 
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Fig.  1 .  Continuous  and  discrete  basis  functions. 

reconstructed  data  for  a=0.59652,N=150 
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Fig.  2.  Transform  data  for  rectangular  pulse 
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Fig.  3.  Reconstructed  signals 


Fig.  4.  The  scattering  structure 
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Fig.  5.  Output  of  the  time  domain  current. 
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